Quantum electrodynamics of a free particle near dispersive dielectric or conducting 
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Quantum electrodynamics near a boundary is investigated by considering the inertial mass shift 
of an electron near a dielectric or conducting surface. We show that in all tractable cases the shift 
can be written in terms of integrals over the TE and TM reflection coefficients associated with the 
surface, in analogy to the Lifshitz formula for the Casimir effect. We discuss the applications and 
potential limitations of this formula, and provide exact results for several models of the surface. 



Quantum electrodynamics is the spectacularly success- 
ful theory of the interaction between charges and electro- 
magnetic fields. The anomalous magnetic moment of the 
electron is perhaps its famous result, finding agreement 
between theory and experiment to an accuracy of one 
part in 10 13 [l|, |2|. However, no conceivable experiment 
can measure a quantity in isolation - there will always be 
apparatus and environment-dependent effects 0-0. One 
type of correction arises from modification of the quan- 
tized electromagnetic field due to material boundaries in 
the vicinity of the system under consideration. Here, we 
consider the effect that the modified quantized field has 
on the self-energy of an electron. While the correspond- 
ing calculations of other quantities such as the magnetic 
moment shift (lfj| have more obvious experimental rele- 
vance, the mass shift calculation turns out to be techni- 
cally simple, and gives exact results even for dispersive 
surfaces. Thus, we present the self energy calculation as 
both an accessible example of our formulation of quan- 
tum field theory near boundaries, and as a clarification 
of a previous result 
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To find the mass shift using the standard formalism of 
quantum field theory, one calculates the self-energy dia- 
gram to the respective order of interest. To 1-loop order 
e 2 = a and in free space, this is a straightforward calcu- 
lation (cf. (l3|)- But, when boundaries are present even 
1-loop calculations of quantum electrodynamics get very 
cumbersome [12| . Since we seek only the change in the 



self energy that is attributable to the surface, we take a 
different, more appropriate approach. While the photon 
propagator receives boundary dependent corrections, the 
electron propagator does not (provided the electron is 
sufficiently far away from the boundary that there is no 
wave-function overlap with the structure of the surface, 
and the interaction between the particle and the surface 
is purely electromagnetic). This means a Feynman di- 
agrammatic approach to the boundary dependent shift 
is not in fact necessary, and it suffices to study a first- 
quantized electron interacting with a second-quantized 
photon field. 

We consider a material filling the space z > 0, de- 
scribed by some dielectric function e(u>). The electron 
sits in vacuum, a distance z = —\z\ away from the ma- 



terial. Via reflection and refraction, the presence of the 
surface affects the electromagnetic field, which interacts 
with the electron via 1 1411 : 



Hmt = p A 

m 



(1) 



where Vi mage is the electrostatic i mag e potential for a 
non-dispersive medium. As shown in |15l | , this interaction 
Hamiltonian remains valid for dispersive media modelled 
as a plasma, via the use of a particular unitary transfor- 
mation. Second-order perturbation theory then gives the 
self energy as 
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for a photon of momentum k and polarization A. We 
make the no-recoil approximation and take the electron's 
final momentum pf to be equal to its initial momentum 
p, which is a reasonable assumption to make as we are 
dealing with a low-energy effect. Writing the quantized 
field A in terms of mode functions fkA and standard pho- 
ton creation and annihilation operators via 



A = E ( f kA fl kA + f kA a k A)> 

all modes 

we can express the self-energy in the form 
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We initially consider two choices of dielectric functions. 
These are a non-dispersive dielectric with e(w) = n 2 , and 
an undamped plasma with dielectric function 
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where uj p is the plasma frequency. The modes fkA can 
be written in terms of plane waves with reflection and 
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transmission coefficients R\ and T\ relevant to incidence 
from either side of the interface. These coefficients are 
given by the standard Fresnel expressions. Writing wave 
vectors as k = (ku,k z ) on the vacuum side and (k»,k d ) 
on the medium side, the Fresnel coefficients that turn 
out to be of importance are the reflection coefficients for 
left-incident radiation 
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k z + k d z * ™ e(uj)k z + k d 

where k d z = ^Je{oj)(k 2 + k 2 ) — k 2 , and we choose 

sgn(fc 2 ) = sgn(fc^). The mode functions in the non- 
dispersive case are: 
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where the modes have been split according to left and right incidence. Barred k vectors correspond to modes that 
have undergone a reflection (reversing the sign of their z component). The modes for the plasma surface are obtained 
from these by the replacements 

/feA, plasma = /fcA.nondisp k 1 ^ e ( u} )) j fkX, plasma = n /fcAaiondisp (" 2 e ( u} )) 

and the addition of a new surfacc-plasmon mode function 
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where we have defined the norming function p(k) = (e 4 — l)/(e 2 \/— 1 — e) and > 0. We split the shift into 
contributions proportional to (p 2 ) and (p 2 ), so that the energy shift is written as 



AE = Y^ ( AE \\ X) + AE 
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Specializing initially to the non-dispersive case, the k z integrals arising in AE^ (or AE^ A> ) can be written in the 
form 
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where V = \/n 2 — ik\t and g\(ku, k z ) is a function an- 
alytic in k z , and which is specific to each polarization. 
The three integrals in the first line represent the left in- 
cident modes, the right-incident travelling modes and the 
right-incident evanescent modes, respectively. The first 
two integrals in the second line are the same as they 
would be in free space, whence we subtract them, since 
we are interested only in the boundary-dependent part of 
the shift. The second two integrals (those in the square 



brackets) can be combined into a single contour integral 
by observing that for k z in the upper half-plane 
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where K > 0, meaning the two integrals can be com- 
bined into 
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The contour C runs over the interval [— oo,0], up the 
imaginary axis to the point ik^y/n 2 — 1/n, back down to 
the real axis and then along [0, — oo]. Taking the same 
approach for the plasma surface, the contour in this case 
is simply [— oo, oo], passing over a branch cut from — u> p to 

I 



uj p . Deforming these contours into the upper half-plane 
results in an additional contribution from a pole at zero 
frequency, and for the plasma case, a contribution from 
a pole in the TM reflection coefficient. We find that this 
latter contribution exactly cancels the surface plasmon 
part of the mode functions (as seen in [ljj 16|). Thus, 
we can write the result for either model in terms of a 
single formula over reflection coefficients, with the choice 
of surface entering only via choice of e(w) in equation ((BJ). 
These integrals are: 
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These formulae arise from residues of a double pole, 
but the dependence on the reflection coefficients is writ- 
ten as explicitly as possible for later convenience. To 
evaluate the results for the two models, one simply inserts 
the relevant dielectric function into the above formulae. 
However, we first consider the case of perfect reflectivity 
- long held to be a good-enough approximation to real 
surfaces. To effect this we take the n 2 — > oo limit of the ^- 



reflection coefficients, resulting in Rf E = — 1, Rtm 
The integrals are trivial, and the result is 



= 1. 
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in agreement with fill ]. Now using e(u) = n 2 in equa- 
tion (|6|), we find 
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The n —> oo limit of this result clearly does not agree 
with the perfect reflector, as discussed in [ll]]. This dis- 
crepancy is one of a family of such issues, we discuss 
the physical meaning and origin of these after presenting 
further results. 

Proceeding, we insert the plasma dielectric function 
(HJ) into (J6j) and carry out the integrals (fT2j) and (fT3|). 
Noting that the k z — > ik\t limits of the plasma reflection 
coefficients coincide with the reflection coefficients for the 
perfect mirror, we find that we can write the plasma shifts 
as corrections to the perfect mirror case via: 
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The integrals Aj_.|| can be evaluated analytically, but 
the resulting expressions arc lengthy combinations of spe- 
cial functions so we will not quote them here. The limit- 



lim A ± ii = 0, 

showing that the plasma and perfect mirror models are 
equivalent at uj p — > oo. 

The surface models discussed so far have all had a de- 
sirable common feature, namely that their corresponding 
wave equations are all Hcrmitian eigenvalue problems (cf. 
[I?} ) ■ A first-principles derivation of the mass shift using 
a dielectric function which does not have this property is 
not possible since one cannot derive the modes. However, 
one can use the Lifshitz theory and write the electromag- 
netic field as a response to fluctuatin g n oise currents in- 
side the material (see, for example, [18|)- The Green's 
function describing this response involves the reflection 
coefficients of the surface so that one necessarily ends up 
with the same formulae (fT2j) and (fT3|). regardless of the 
specific choice of dielectric function. Proceeding along 
these lines, we use our formulae to investigate a model 
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in which a restoring force is introduced [19( , summarized 
by the introduction of a new characteristic frequency lot 
into the dielectric function: 



Surprisingly, the integrals are here much simpler than 
for the plasma case, the result is: 



e( W ) 



(18) 



AE] 



Disp 



e 2 ^ 



167rm 2 (lo 2 + 2lo 2 Y 



1 

2z* 



(pi 



- + -(2^ 



We can compare this result to the non-dispersive case 
by using the static susceptibility x(0) 



X(0) = e(0) - 1 



n 2 — 1 non-dispersive 
uip I lo t dispersive 



This gives, for example, the perpendicular component 



as 
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It is seen that, for large x(0)j the dispersive dielectric 
is in agreement with the non-dispersive dielectric with 
large n. Figure [1] shows the energy shift as dependent on 
the static susceptibility for lotz = 0.2. The plasma result 
cannot be shown in Fig. ([T]) since its static susceptibility 
is infinite; however it can be compared to the dispersive 
shift via a plot against the dimensionless parameter lo p z, 
as shown in Fig. ([2]) . In both plots we show energy shifts 
in units of the perfect reflector shift. 
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FIG. 1. [Colour online] Energy shift vs static susceptibility 
for various models in units of the perfect reflector shift. The 
position and height of the peak may be tuned by varying 
lotz, and the factor-of-two disagreement between the perfect 
reflector and dielectric cases is shown. 




FIG. 2. [Colour online] Energy shift in units of the perfect re- 
flector shift vs the dimensionless parameter u) p z for the plasma 
and dispersive dielectric models for various lotz. The values 
of ujtz shown are 0.2 (solid line), 0.4 (dashed) and 0.6 (dot- 
ted). 



It is not hard to show that the peak emerges only when 
lutz < 1/y/E. It moves towards x(0) = 2 for decreasing 
lotz, and its height scales as (lotz)~ 2 . Thus, for small 
values of this lotz we see that the shift can be made con- 
siderably larger than in the previously considered perfect 
reflector model. 

Now we turn our attention to the origins of the dis- 
crepancies between the models. Mathematically, the dis- 
agreements arise because of non-commutation of limits 
in the reflection coefficients (or their derivatives), namely 
k z — > ikit (required for finding the residue at this point) 
and whatever limit one has to take to get from dielectric 
function to another. For example, the lot limit of 
the result for the dispersive dielectric should take us to 
the plasma result, but it does not. This is because the 
lot — > and k z — > ik\\ limits of the derivative of the TM 
reflection coefficient do not commute. A similar problem 
causes the perfect reflector and non-dispersive results to 
disagree in the limit n — > oo. 

Physically, the differences between models that dis- 
agree with each other are down to a number of reasons. 
One of them is the exclusion of part of the photon phase- 
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space, namely the evanescent modes. Previous workers 
have shown that exclusion of evanescent modes is not an 



k x is: 



adequate approximation to reality J12J, we confirm this 
conclusion in the context of a dispersive medium. The 
other main reason has to do with the different response 
of conductors and dielectrics to electric fields at low fre- 
quencies: e(w) has a pole at lo = for a conductor but 
not for a dielectric. The discrepancies between the results 
for the mass shift show that one has to decide whether 
the material at hand should be modelled as a metal (no 
restoring force for the charge carriers) or as a dielectric 
(with a restoring force parametrized by u>t), since these 
two classes of model for the surface are not obtainable 
as limiting cases of one another, reflecting the different 
nature of the electromagnetic response of conductors and 
dielectrics. 

To find conductors and dielectrics giving rise to differ- 
ent results on account of their different response to elec- 
tromagnetic fields is of course not at all surprising. This 
is, however, in contrast to what one might have expected 
from the closely related Casimir-Polder energy shift in 
an atom close to a conducting or dielectric boundary: 
in both retarded and non-retarded regimes the Casimir- 
Polder shift of an atom in front of a dielectric [2(| repro- 
duces the original result for an atom close to a perfect re- 
flector [2l[ in the limit of infinite dielectric constant, and 
so does the level shift for an atom near a plasma surface 
22j in the limit of infinite plasma frequency, ui p — > oo. 
The crucial difference between an atom and a free particle 
in this context is that the excitation spectrum of a bound 
electron has a gap at low frequencies corresponding to the 
nearest energy level whereas a free particle admit exci- 
tations of arbitrarily low frequency. As a consequence, 
the low-frequency behaviour of the electromagnetic re- 
sponse of the material, in particular the pole at ui = 
in the dielectric function of a conductor, play a decisive 
role for the mass shift of a free particle, but not for the 
Casimir-Polder shift of an atom. 

The decisive importance of the pole at lo = in the 
dielectric function of a conductor is made obvious by the 
fact that the energy shifts (|16p and (| 1 T[) do not vanish 



in the limit 



0, despite e(w) reducing to the vac- 



uum value of 1 in that case. The limit ui p — > is non- 
analytic because the choice of a dielectric function of the 
form ([5]) necessarily describes freely moving charge carri- 
ers at w = 0, which is obviously not true for vacuum with 
e=l. Mathematically speaking, eq. ([5]) is ill-defined if 
both uj — > and ui p — > 0; in line with the physical inter- 
pretation, the fact that e(ui) has a pole at us = is more 
important than the strength of this pole. 



A natural next step in the investigation of the effect 
for realistic materials would be to include a damping pa- 
rameter, meaning that the dielectric function in terms of 
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This introduces the additional complication that the re- 
flection coefficient has branch points at k z = iifcii , caus- 
ing the formulae (fT2"| and ([TS]) to become ambiguous. 
Our method is reliant on the fact that whatever happens 
below iku in the upper half of the k z plane does not pre- 
clude a contour deformation that allows us to simply pick 
up the residue at k z = iku . This inability to deform a con- 
tour once damping is introduced has been investigated in 
detail in [23| , where it is concluded that standard deriva- 
tions of the Lifshitz formula with damping must contain 
an unspecified inconsistency (a statement reflected by ex- 
perimental results [24]]). In our approach, the failure of 
the damped reflection coefficient to be single-valued in 
the neighborhood of k z = ik\\ could be a manifestation 
of these problems with damping detailed by others. 

Our results are intimately related to the shift in the 
cyclotron frequency of an electron near a surface. If the 
external magnetic field is directed perpendicular to the 
surface, the calculations coincide, so that a measurement 
of the cyclotron frequency is in effect a measurement of 
the mass shift. Precision g — 2 experiments often rely on 
accurate measurement of the cyclotron frequency [3, 0] , 
but a much more elaborate calculation would have to be 
undertaken in order to precisely enumerate the effects 
relevant to a specific apparatus. Finally we note that 
for magnetic fields directed parallel to the surface, the 
additional electrostatic interaction skews the orbit, and 
much more so than the mass anisotropy Thus, a 
measurement of the cyclotron frequency in a parallel field 
does not deliver the mass shift. 

In summary, we have calculated the self-energy of an 
electron near dispersive surfaces of various kinds, which 
had previously been considered only in idealized mod- 
els. We have shown that is it crucial that one decides at 
the start of a calculation whether the material should be 
modelled as a metal or a dielectric, since the results for 
the two classes of material are not obtainable as limiting 
cases of one another. 
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